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GENERATING AND DISPLAYING SPATIALLY OFFSET SUB- 
FRAMES ON DIFFERENT TYPES OF GRIDS 

Cross-Reference to Related Applications 

5 This application is related to U.S. Patent Application Serial No. 

10/213,555, filed on August 7, 2002, entitled IMAGE DISPLAY SYSTEM 
AND METHOD; U.S. Patent Application Serial No. 10/242,195, filed on 
September 1 1, 2002, entitled IMAGE DISPLAY SYSTEM AND METHOD; 
U.S. Patent Application Serial No. 10/242,545, filed on September 11, 2002, 

10 entitled IMAGE DISPLAY SYSTEM AND METHOD; U.S. Patent Application 
Serial No. 10/631,681, filed July 31, 2003, entitled GENERATING AND 
DISPLAYING SPATIALLY OFFSET SUB-FRAMES; U.S. Patent Application 
Serial No. 10/632,042, filed July 31, 2003, entitled GENERATING AND 
DISPLAYING SPATIALLY OFFSET SUB-FRAMES; U.S. Patent Application 

15 Serial No. 10/672,845, filed September 26, 2003, entitled GENERATING AND 
DISPLAYING SPATIALLY OFFSET SUB-FRAMES; U.S. Patent Application 
Serial No. 10/672,544, filed September 26, 2003, entitled GENERATING AND 
DISPLAYING SPATIALLY OFFSET SUB-FRAMES; U.S. Patent Application 
Serial No. , attorney docket number 200314247-1, filed on the 

20 same date as the present application, and entitled GENERATING AND 

DISPLAYING SPATIALLY OFFSET SUB-FRAMES ON A DIAMOND 

GRID; and U.S. Patent Application Serial No. , attorney docket 

number 10015701-1, filed on the same date as the present application, and 
entitled IMAGE DISPLAY SYSTEM AND METHOD. Each of the above U.S. 

25 Patent Applications is assigned to the assignee of the present invention, and is 
hereby incorporated by reference herein. 

The Field of the Invention 

The present invention generally relates to display systems, and more 
30 particularly to generating and displaying spatially offset sub-frames on different 
types of grids. 
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Background of the Invention 

A conventional system or device for displaying an image, such as a 
display, projector, or other imaging system, produces a displayed image by 
addressing an array of individual picture elements or pixels arranged in a pattern, 
5 such as in horizontal rows and vertical columns, a diamond grid, or other pattern. 
A resolution of the displayed image for a pixel pattern with horizontal rows and 
vertical columns is defined as the number of horizontal rows and vertical 
columns of individual pixels forming the displayed image. The resolution of the 
displayed image is affected by a resolution of the display device itself as well as 

10 a resolution of the image data processed by the display device and used to 
produce the displayed image. 

Typically, to increase a resolution of the displayed image, the resolution 
of the display device as well as the resolution of the image data used to produce 
the displayed image must be increased. Increasing a resolution of the display 

15 device, however, increases a cost and complexity of the display device. In 

addition, higher resolution image data may not be available or may be difficult to 
generate. 

Summary of the Invention 

20 One form of the present invention provides a method of displaying an 

image with a display device. The method includes receiving image data for the 
image on a first type of grid. The method includes generating a first sub-frame 
and a second sub-frame corresponding to the image data, the first and the second 
sub-frames each generated on a second type of grid that is different than the first 

25 type of grid. The method includes alternating between displaying the first sub- 
frame in a first position and displaying the second sub-frame in a second position 
spatially offset from the first position. 

Brief Description of the Drawings 

30 Figure 1 is a block diagram illustrating an image display system 

according to one embodiment of the present invention. 
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Figures 2A-2C are schematic diagrams illustrating the display of two 
sub-frames according to one embodiment of the present invention. 

Figures 3A-3E are schematic diagrams illustrating the display of four 
sub-frames according to one embodiment of the present invention. 
5 Figures 4A-4E are schematic diagrams illustrating the display of a pixel 

with an image display system according to one embodiment of the present 
invention. 

Figure 5 is a diagram illustrating the generation of low resolution sub- 
frames from an original high resolution image using a nearest neighbor 
10 algorithm according to one embodiment of the present invention. 

Figure 6 is a diagram illustrating the generation of low resolution sub- 
frames from an original high resolution image using a bilinear algorithm 
according to one embodiment of the present invention. 

Figure 7 is a block diagram illustrating a system for generating a 
15 simulated high resolution image according to one embodiment of the present 
invention. 

Figure 8 is a block diagram illustrating a system for generating a 
simulated high resolution image for two-position processing based on separable 
upsampling according to one embodiment of the present invention. 
20 Figure 9 is a block diagram illustrating a system for generating a 

simulated high resolution image for two-position processing based on non- 
separable upsampling according to one embodiment of the present invention. 

Figure 10 is a block diagram illustrating a system for generating a 
simulated high resolution image for four-position processing according to one 
25 embodiment of the present invention. 

Figure 1 1 is a block diagram illustrating the comparison of a simulated 
high resolution image and a desired high resolution image according to one 
embodiment of the present invention. 

Figure 12 is a diagram illustrating the effect in the frequency domain of 
30 the upsampling of a sub-frame according to one embodiment of the present 
invention. 



3 



PATENT 
200314885-1 



Figure 13 is a diagram illustrating the effect in the frequency domain of 
the shifting of an upsampled sub-frame according to one embodiment of the 
present invention. 

Figure 14 is a diagram illustrating regions of influence for pixels in an 
5 upsampled image according to one embodiment of the present invention. 

Figure 15 is a diagram illustrating the generation of an initial simulated 
high resolution image based on an adaptive multi-pass algorithm according to 
one embodiment of the present invention. 

Figure 16 is a diagram illustrating the generation of correction data based 
10 on an adaptive multi-pass algorithm according to one embodiment of the present 
invention. 

Figure 17 is a diagram illustrating the generation of updated sub-frames 
based on an adaptive multi-pass algorithm according to one embodiment of the 
present invention. 

15 Figure 18 is a diagram illustrating the generation of correction data based 

on an adaptive multi-pass algorithm according to another embodiment of the 
present invention. 

Figure 19A is a diagram illustrating rectangular-shaped pixels on a 
rectangular grid according to one embodiment of the present invention. 

20 Figure 19B is a diagram illustrating diamond-shaped pixels on a diamond 

grid according to one embodiment of the present invention. 

Figure 20 is a diagram illustrating the display of two sub-frames with 
rectangular-shaped pixels on rectangular grids and a diagonal offset between 
sub-frames according to one embodiment of the present invention. 

25 Figure 21 is a diagram illustrating the display of two sub-frames with 

diamond-shaped pixels on diamond grids and a horizontal offset between sub- 
frames according to one embodiment of the present invention. 

Figure 22 is a block diagram illustrating a system for generating a 
simulated high resolution image on a rectangular grid based on sub-frames on 

30 diamond grids according to one embodiment of the present invention. 



4 



o 



10 



PATENT 
200314885-1 



Figure 23 is a block diagram illustrating a system for generating a 
simulated high resolution image on a rectangular grid based on an image on a 
diamond grid according to one embodiment of the present invention. 

Figure 24 is a diagram illustrating the display of two sub-frames with 
rectangular-shaped pixels on rectangular grids and a diagonal offset between 
sub-frames according to one embodiment of the present invention. 

Figure 25 is a block diagram illustrating a system for generating a 
simulated high resolution image on a diamond grid based on sub-frames on 
rectangular grids according to one embodiment of the present invention. 



Description of the Preferred Embodiments 

In the following detailed description of the preferred embodiments, 
reference is made to the accompanying drawings, which form a part hereof, and 
in which is shown by way of illustration specific embodiments in which the 
15 invention may be practiced. It is to be understood that other embodiments may 
be utilized and structural or logical changes may be made without departing 
from the scope of the present invention. The following detailed description, 
therefore, is not to be taken in a limiting sense, and the scope of the present 
invention is defined by the appended claims. 

20 

L Spatial and Temporal Shifting of Sub-frames 

Some display systems, such as some digital light projectors, may not 
have sufficient resolution to display some high resolution images. Such systems 
can be configured to give the appearance to the human eye of higher resolution 

25 images by displaying spatially and temporally shifted lower resolution images. 
The lower resolution images are referred to as sub-frames. A problem of sub- 
frame generation, which is addressed by embodiments of the present invention, 
is to determine appropriate values for the sub-frames so that the displayed sub- 
frames are close in appearance to how the high-resolution image from which the 

30 sub-frames were derived would appear if directly displayed. 

One embodiment of a display system that provides the appearance of 
enhanced resolution through temporal and spatial shifting of sub-frames is 
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described in the above-cited U.S. patent applications, and is summarized below 
with reference to Figures 1-4E. 

Figure 1 is a block diagram illustrating an image display system 10 
according to one embodiment of the present invention. Image display system 10 
5 facilitates processing of an image 12 to create a displayed image 14. Image 12 is 
defined to include any pictorial, graphical, or textural characters, symbols, 
illustrations, or other representation of information. Image 12 is represented, for 
example, by image data 16. Image data 16 includes individual picture elements 
or pixels of image 12. While one image is illustrated and described as being 

10 processed by image display system 10, it is understood that a plurality or series 
of images may be processed and displayed by image display system 10. 

In one embodiment, image display system 10 includes a frame rate 
conversion unit 20 and an image frame buffer 22, an image processing unit 24, 
and a display device 26. As described below, frame rate conversion unit 20 and 

15 image frame buffer 22 receive and buffer image data 16 for image 12 to create 
an image frame 28 for image 12. Image processing unit 24 processes image 
frame 28 to define one or more image sub-frames 30 for image frame 28, and 
display device 26 temporally and spatially displays image sub-frames 30 to 
produce displayed image 14. 

20 Image display system 10, including frame rate conversion unit 20 and 

image processing unit 24, includes hardware, software, firmware, or a 
combination of these. In one embodiment, one or more components of image 
display system 10, including frame rate conversion unit 20 and image processing 
unit 24, are included in a computer, computer server, or other microprocessor- 

25 based system capable of performing a sequence of logic operations. In addition, 
processing can be distributed throughout the system with individual portions 
being implemented in separate system components. 

Image data 16 may include digital image data 161 or analog image data 
162. To process analog image data 162, image display system 10 includes an 

30 analog-to-digital (A/D) converter 32. As such, A/D converter 32 converts 

analog image data 162 to digital form for subsequent processing. Thus, image 
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display system 10 may receive and process digital image data 161 or analog 
image data 162 for image 12. 

Frame rate conversion unit 20 receives image data 16 for image 12 and 
buffers or stores image data 16 in image frame buffer 22. More specifically, 
5 frame rate conversion unit 20 receives image data 16 representing individual 

lines or fields of image 12 and buffers image data 16 in image frame buffer 22 to 
create image frame 28 for image 12. Image frame buffer 22 buffers image data 
16 by receiving and storing all of the image data for image frame 28, and frame 
rate conversion unit 20 creates image frame 28 by subsequently retrieving or 

10 extracting all of the image data for image frame 28 from image frame buffer 22. 
As such, image frame 28 is defined to include a plurality of individual lines or 
fields of image data 16 representing an entirety of image 12. In one 
embodiment, image frame 28 includes a plurality of columns and a plurality of 
rows of individual pixels on a rectangular grid representing image 12. In another 

15 embodiment, image frame 28 includes a plurality of pixels on a diamond grid 
representing image 12. 

Frame rate conversion unit 20 and image frame buffer 22 can receive and 
process image data 16 as progressive image data or interlaced image data. With 
progressive image data, frame rate conversion unit 20 and image frame buffer 22 

20 receive and store sequential fields of image data 16 for image 12. Thus, frame 
rate conversion unit 20 creates image frame 28 by retrieving the sequential fields 
of image data 16 for image 12. With interlaced image data, frame rate 
conversion unit 20 and image frame buffer 22 receive and store odd fields and 
even fields of image data 16 for image 12. For example, all of the odd fields of 

25 image data 16 are received and stored and all of the even fields of image data 16 
are received and stored. As such, frame rate conversion unit 20 de-interlaces 
image data 16 and creates image frame 28 by retrieving the odd and even fields 
of image data 16 for image 12. 

Image frame buffer 22 includes memory for storing image data 16 for 

30 one or more image frames 28 of respective images 12. Thus, image frame buffer 
22 constitutes a database of one or more image frames 28. Examples of image 
frame buffer 22 include non-volatile memory (e.g., a hard disk drive or other 
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persistent storage device) and may include volatile memory (e.g., random access 
memory (RAM)). 

By receiving image data 16 at frame rate conversion unit 20 and 
buffering image data 16 with image frame buffer 22, input timing of image data 
5 16 can be decoupled from a timing requirement of display device 26. More 
specifically, since image data 16 for image frame 28 is received and stored by 
image frame buffer 22, image data 16 can be received as input at any rate. As 
such, the frame rate of image frame 28 can be converted to the timing 
requirement of display device 26. Thus, image data 16 for image frame 28 can 

10 be extracted from image frame buffer 22 at a frame rate of display device 26. 

In one embodiment, image processing unit 24 includes a resolution 
adjustment unit 34 and a sub-frame generation unit 36. As described below, 
resolution adjustment unit 34 receives image data 16 for image frame 28 and 
adjusts a resolution of image data 16 for display on display device 26, and sub- 

15 frame generation unit 36 generates a plurality of image sub-frames 30 for image 
frame 28. More specifically, image processing unit 24 receives image data 16 
for image frame 28 at an original resolution and processes image data 16 to 
increase, decrease, or leave unaltered the resolution of image data 16. 
Accordingly, with image processing unit 24, image display system 10 can 

20 receive and display image data 16 of varying resolutions. 

Sub-frame generation unit 36 receives and processes image data 16 for 
image frame 28 to define a plurality of image sub-frames 30 for image frame 28. 
If resolution adjustment unit 34 has adjusted the resolution of image data 16, 
sub-frame generation unit 36 receives image data 16 at the adjusted resolution. 

25 The adjusted resolution of image data 16 may be increased, decreased, or the 

same as the original resolution of image data 16 for image frame 28. Sub-frame 
generation unit 36 generates image sub-frames 30 with a resolution which 
matches the resolution of display device 26. Image sub-frames 30 are each of an 
area equal to image frame 28. In one embodiment, sub-frames 30 each include a 

30 plurality of columns and a plurality of rows of individual pixels on a rectangular 
grid representing a subset of image data 16 of image 12. In another 
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embodiment, sub-frames 30 each include a plurality of pixels arranged on a 
diamond grid. 

Image sub-frames 30 are spatially offset from each other when displayed. 
In one embodiment, image sub-frames 30 are offset from each other by a vertical 
5 distance and a horizontal distance, as described below. 

Display device 26 receives image sub-frames 30 from image processing 
unit 24 and sequentially displays image sub-frames 30 to create displayed image 
14. More specifically, as image sub-frames 30 are spatially offset from each 
other, display device 26 displays image sub-frames 30 in different positions 

10 according to the spatial offset of image sub-frames 30, as described below. As 
such, display device 26 alternates between displaying image sub-frames 30 for 
image frame 28 to create displayed image 14. Accordingly, display device 26 
displays an entire sub-frame 30 for image frame 28 at one time. 

In one embodiment, display device 26 performs one cycle of displaying 

15 image sub-frames 30 for each image frame 28. Display device 26 displays 

image sub-frames 30 so as to be spatially and temporally offset from each other. 
In one embodiment, display device 26 optically steers image sub-frames 30 to 
create displayed image 14. As such, individual pixels of display device 26 are 
addressed to multiple locations. 

20 In one embodiment, display device 26 includes an image shifter 38. 

Image shifter 38 spatially alters or offsets the position of image sub-frames 30 as 
displayed by display device 26. More specifically, image shifter 38 varies the 
position of display of image sub-frames 30, as described below, to produce 
displayed image 14. 

25 In one embodiment, display device 26 includes a light modulator for 

modulation of incident light. The light modulator includes, for example, a 
plurality of micro-mirror devices arranged to form an array of micro-mirror 
devices. As such, each micro-mirror device constitutes one cell or pixel of 
display device 26. Display device 26 may form part of a display, projector, or 

30 other imaging system. 

In one embodiment, image display system 10 includes a timing generator 
40. Timing generator 40 communicates, for example, with frame rate 
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conversion unit 20, image processing unit 24, including resolution adjustment 
unit 34 and sub-frame generation unit 36, and display device 26, including image 
shifter 38. As such, timing generator 40 synchronizes buffering and conversion 
of image data 16 to create image frame 28, processing of image frame 28 to 
5 adjust the resolution of image data 16 and generate image sub-frames 30, and 
positioning and displaying of image sub-frames 30 to produce displayed image 
14. Accordingly, timing generator 40 controls timing of image display system 
10 such that entire sub-frames of image 12 are temporally and spatially 
displayed by display device 26 as displayed image 14. 

10 In one embodiment, as illustrated in Figures 2 A and 2B, image 

processing unit 24 defines two image sub-frames 30 for image frame 28. More 
specifically, image processing unit 24 defines a first sub-frame 301 and a second 
sub-frame 302 for image frame 28. As such, first sub-frame 301 and second 
sub-frame 302 each include a plurality of columns and a plurality of rows of 

15 individual pixels 18 of image data 16. Thus, first sub-frame 301 and second sub- 
frame 302 each constitute an image data array or pixel matrix of a subset of 
image data 16. 

In one embodiment, as illustrated in Figure 2B, second sub-frame 302 is 
offset from first sub-frame 301 by a vertical distance 50 and a horizontal 
20 distance 52. As such, second sub-frame 302 is spatially offset from first sub- 
frame 301 by a predetermined distance. In one illustrative embodiment, vertical 
distance 50 and horizontal distance 52 are each approximately one-half of one 
pixel. 

As illustrated in Figure 2C, display device 26 alternates between 
25 displaying first sub-frame 301 in a first position and displaying second sub- 
frame 302 in a second position spatially offset from the first position. More 
specifically, display device 26 shifts display of second sub-frame 302 relative to 
display of first sub-frame 301 by vertical distance 50 and horizontal distance 52. 
As such, pixels of first sub-frame 301 overlap pixels of second sub-frame 302. 
30 In one embodiment, display device 26 performs one cycle of displaying first sub- 
frame 301 in the first position and displaying second sub-frame 302 in the 
second position for image frame 28. Thus, second sub-frame 302 is spatially and 

10 
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temporally displayed relative to first sub-frame 301. The display of two 
temporally and spatially shifted sub-frames in this manner is referred to herein as 
two-position processing. 

In another embodiment, as illustrated in Figures 3A-3D, image 
5 processing unit 24 defines four image sub-frames 30 for image frame 28. More 
specifically, image processing unit 24 defines a first sub-frame 301, a second 
sub-frame 302, a third sub-frame 303, and a fourth sub-frame 304 for image 
frame 28. As such, first sub-frame 301, second sub-frame 302, third sub-frame 
303, and fourth sub-frame 304 each include a plurality of columns and a 

10 plurality of rows of individual pixels 18 of image data 16. 

In one embodiment, as illustrated in Figures 3B-3D, second sub-frame 
302 is offset from first sub-frame 301 by a vertical distance 50 and a horizontal 
distance 52, third sub-frame 303 is offset from first sub-frame 301 by a 
horizontal distance 54, and fourth sub-frame 304 is offset from first sub-frame 

15 301 by a vertical distance 56. As such, second sub-frame 302, third sub-frame 
303, and fourth sub-frame 304 are each spatially offset from each other and 
spatially offset from first sub-frame 301 by a predetermined distance. In one 
illustrative embodiment, vertical distance 50, horizontal distance 52, horizontal 
distance 54, and vertical distance 56 are each approximately one-half of one 

20 pixel. 

As illustrated schematically in Figure 3E, display device 26 alternates 
between displaying first sub-frame 301 in a first position Pi, displaying second 
sub-frame 302 in a second position P 2 spatially offset from the first position, 
displaying third sub-frame 303 in a third position P 3 spatially offset from the 

25 first position, and displaying fourth sub-frame 304 in a fourth position P 4 
spatially offset from the first position. More specifically, display device 26 
shifts display of second sub-frame 302, third sub-frame 303, and fourth sub- 
frame 304 relative to first sub-frame 301 by the respective predetermined 
distance. As such, pixels of first sub-frame 301, second sub-frame 302, third 

30 sub-frame 303, and fourth sub-frame 304 overlap each other. 

In one embodiment, display device 26 performs one cycle of displaying 
first sub-frame 301 in the first position, displaying second sub-frame 302 in the 

11 
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second position, displaying third sub-frame 303 in the third position, and 
displaying fourth sub-frame 304 in the fourth position for image frame 28. 
Thus, second sub-frame 302, third sub-frame 303, and fourth sub-frame 304 are 
spatially and temporally displayed relative to each other and relative to first sub- 
5 frame 301. The display of four temporally and spatially shifted sub-frames in 
this manner is referred to herein as four-position processing. 

Figures 4A-4E illustrate one embodiment of completing one cycle of 
displaying a pixel 181 from first sub-frame 301 in the first position, displaying a 
pixel 182 from second sub-frame 302 in the second position, displaying a pixel 

10 183 from third sub-frame 303 in the third position, and displaying a pixel 184 
from fourth sub-frame 304 in the fourth position. More specifically, Figure 4A 
illustrates display of pixel 181 from first sub-frame 301 in the first position, 
Figure 4B illustrates display of pixel 182 from second sub-frame 302 in the 
second position (with the first position being illustrated by dashed lines), Figure 

15 4C illustrates display of pixel 183 from third sub-frame 303 in the third position 
(with the first position and the second position being illustrated by dashed lines), 
Figure 4D illustrates display of pixel 184 from fourth sub-frame 304 in the 
fourth position (with the first position, the second position, and the third position 
being illustrated by dashed lines), and Figure 4E illustrates display of pixel 181 

20 from first sub-frame 301 in the first position (with the second position, the third 
position, and the fourth position being illustrated by dashed lines). 

Sub-frame generation unit 36 (Figure 1) generates sub-frames 30 based 
on image data in image frame 28. It will be understood by a person of ordinary 
skill in the art that functions performed by sub-frame generation unit 36 may be 

25 implemented in hardware, software, firmware, or any combination thereof. The 
implementation may be via a microprocessor, programmable logic device, or 
state machine. Components of the present invention may reside in software on 
one or more computer-readable mediums. The term computer-readable medium 
as used herein is defined to include any kind of memory, volatile or non-volatile, 

30 such as floppy disks, hard disks, CD-ROMs, flash memory, read-only memory 
(ROM), and random access memory. 
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In one form of the invention, sub-frames 30 have a lower resolution than 
image frame 28. Thus, sub-frames 30 are also referred to herein as low 
resolution images 30, and image frame 28 is also referred to herein as a high 
resolution image 28. It will be understood by persons of ordinary skill in the art 
5 that the terms low resolution and high resolution are used herein in a 

comparative fashion, and are not limited to any particular minimum or maximum 
number of pixels. In one embodiment, sub-frame generation unit 36 is 
configured to generate sub-frames 30 based on one of five algorithms. These 
five algorithms are referred to herein as the following: (1) nearest neighbor; (2) 

10 bilinear; (3) spatial domain; (4) frequency domain; and (5) adaptive multi-pass. 

The nearest neighbor algorithm and the bilinear algorithm according to 
one form of the invention generate sub-frames 30 by combining pixels from a 
high resolution image 28. The spatial domain algorithm and the frequency 
domain algorithm according to one form of the invention generate sub-frames 30 

15 based on the minimization of a global error metric that represents a difference 
between a simulated high resolution image and a desired high resolution image 
28. The adaptive multi-pass algorithm according to one form of the invention 
generates sub-frames 30 based on the minimization of a local error metric. In 
one embodiment, sub-frame generation unit 36 includes memory for storing a 

20 relationship between sub-frame values and high resolution image values, 

wherein the relationship is based on minimization of an error metric between the 
high resolution image values and a simulated high resolution image that is a 
function of the sub-frame values. Embodiments of each of these five algorithms 
are described below with reference to Figures 5-18. 

25 

II. Nearest Neighbor 

Figure 5 is a diagram illustrating the generation of low resolution sub- 
frames 30A and 30B from an original high resolution image 28 using a nearest 
neighbor algorithm according to one embodiment of the present invention. In 
30 the illustrated embodiment, high resolution image 28 includes four columns and 
four rows of pixels, for a total of sixteen pixels HI -HI 6. In one embodiment of 
the nearest neighbor algorithm, a first sub-frame 30A is generated by taking 

13 
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every other pixel in a first row of the high resolution image 28, skipping the 
second row of the high resolution image 28, taking every other pixel in the third 
row of the high resolution image 28, and repeating this process throughout the 
high resolution image 28. Thus, as shown in Figure 5, the first row of sub-frame 
5 30A includes pixels HI and H3, and the second row of sub-frame 30A includes 
pixels H9 and HI 1. In one form of the invention, a second sub-frame 30B is 
generated in the same manner as the first sub-frame 30A, but the process begins 
at a pixel H6 that is shifted down one row and over one column from the first 
pixel HI. Thus, as shown in Figure 5, the first row of sub-frame 30B includes 
10 pixels H6 and H8, and the second row of sub-frame 30B includes pixels H14 and 
H16. 

In one embodiment, the nearest neighbor algorithm is implemented with 
a 2x2 filter with three filter coefficients of "0" and a fourth filter coefficient of 
"1" to generate a weighted sum of the pixel values from the high resolution 
15 image. Displaying sub-frames 30A and 30B using two-position processing as 
described above gives the appearance of a higher resolution image. The nearest 
neighbor algorithm is also applicable to four-position processing, and is not 
limited to images having the number of pixels shown in Figure 5. 

20 HI. Bilinear 

Figure 6 is a diagram illustrating the generation of low resolution sub- 
frames 30C and 30D from an original high resolution image 28 using a bilinear 
algorithm according to one embodiment of the present invention. In the 
illustrated embodiment, high resolution image 28 includes four columns and four 

25 rows of pixels, for a total of sixteen pixels HI -HI 6. Sub-frame 30C includes 
two columns and two rows of pixels, for a total of four pixels L1-L4. And sub- 
frame 30D includes two columns and two rows of pixels, for a total of four 
pixels L5-L8. 

In one embodiment, the values for pixels L1-L8 in sub-frames 30C and 
30 30D are generated from the pixel values HI -HI 6 of image 28 based on the 
following Equations I-VIII: 
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Equation I 

LI =(4H1 +2H2 +2H5)/8 
Equation II 

L2 = (4H3 + 2H4 +2H7) / 8 
5 Equation III 

L3 = (4H9 + 2H10 +2H13) / 8 
Equation IV 

L4 = (4H11 +2H12+2H15)/8 
Equation V 

10 L5 = (4H6 + 2H2 +2H5) / 8 

Equation VI 

L6 = (4H8 + 2H4 +2H7) / 8 
Equation VII 

L7 = (4H14 + 2H10 +2H13) / 8 
15 Equation VIII 

L8 = (4H16 + 2H12 +2H15) / 8 
As can be seen from the above Equations I- VIE, the values of the pixels 
L1-L4 in sub-frame 30C are influenced the most by the values of pixels HI, H3, 
H9, and HI 1, respectively, due to the multiplication by four. But the values for 
20 the pixels L1-L4 in sub-frame 30C are also influenced by the values of diagonal 
neighbors of pixels HI, H3, H9, and HI 1. Similarly, the values of the pixels L5- 
L8 in sub-frame 30D are influenced the most by the values of pixels H6, H8, 
H14, and HI 6, respectively, due to the multiplication by four. But the values for 
the pixels L5-L8 in sub-frame 30D are also influenced by the values of diagonal 
25 neighbors of pixels H6, H8, H14, and H16. 

In one embodiment, the bilinear algorithm is implemented with a 2x2 
filter with one filter coefficient of "0" and three filter coefficients having a non- 
zero value (e.g., 4, 2, and 2) to generate a weighted sum of the pixel values from 
the high resolution image. In another embodiment, other values are used for the 
30 filter coefficients. Displaying sub-frames 30C and 30D using two-position 

processing as described above gives the appearance of a higher resolution image. 
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The bilinear algorithm is also applicable to four-position processing, and is not 
limited to images having the number of pixels shown in Figure 6. 

In one form of the nearest neighbor and bilinear algorithms, sub-frames 
30 are generated based on a linear combination of pixel values from an original 
5 high resolution image as described above. In another embodiment, sub-frames 
30 are generated based on a non-linear combination of pixel values from an 
original high resolution image. For example, if the original high resolution 
image is gamma-corrected, appropriate non-linear combinations are used in one 
embodiment to undo the effect of the gamma curve. 

10 

IV. Systems for Generating Simulated High Resolution Images 

Figures 7-10 illustrate systems for generating simulated high resolution 
images. Based on these systems, spatial domain, frequency domain, and 
adaptive multi-pass algorithms for generating sub-frames are developed, as 

15 described in further detail below. 

Figure 7 is a block diagram illustrating a system 400 for generating a 
simulated high resolution image 412 from two 4x4 pixel low resolution sub- 
frames 30E according to one embodiment of the present invention. System 400 
includes upsampling stage 402, shifting stage 404, convolution stage 406, and 

20 summation stage 410. Sub-frames 30E are upsampled by upsampling stage 402 
based on a sampling matrix, M, thereby generating upsampled images. The 
upsampled images are shifted by shifting stage 404 based on a spatial shifting 
matrix, S, thereby generating shifted upsampled images. The shifted upsampled 
images are convolved with an interpolating filter at convolution stage 406, 

25 thereby generating blocked images 408. In the illustrated embodiment, the 
interpolating filter is a 2x2 filter with filter coefficients of "1", and with the 
center of the convolution being the upper left position in the 2x2 matrix. The 
interpolating filter simulates the superposition of low resolution sub-frames on a 
high resolution grid. The low resolution sub-frame pixel data is expanded so 

30 that the sub-frames can be represented on a high resolution grid. The 

interpolating filter fills in the missing pixel data produced by upsampling. The 



16 



PATENT 
200314885-1 

blocked images 408 are weighted and summed by summation block 410 to 
generate the 8x8 pixel simulated high resolution image 412. 

Figure 8 is a block diagram illustrating a system 500 for generating a 
simulated high resolution image 512 for two-position processing based on 
5 separable upsampling of two 4x4 pixel low resolution sub-frames 30F and 30G 
according to one embodiment of the present invention. System 500 includes 
upsampling stages 502 and 514, shifting stage 518, convolution stages 506 and 
522, summation stage 508, and multiplication stage 510. Sub-frame 30F is 
upsampled by a factor of two by upsampling stage 502, thereby generating an 

10 8x8 pixel upsampled image 504. The dark pixels in upsampled image 504 
represent the sixteen pixels from sub-frame 30F, and the light pixels in 
upsampled image 504 represent zero values. Sub-frame 30G is upsampled by a 
factor of two by upsampling stage 514, thereby generating an 8x8 pixel 
upsampled image 516. The dark pixels in upsampled image 516 represent the 

15 sixteen pixels from sub-frame 30G, and the light pixels in upsampled image 516 
represent zero values. In one embodiment, upsampling stages 502 and 514 
upsample sub-frames 30F and 30G, respectively, using a diagonal sampling 
matrix. 

The upsampled image 516 is shifted by shifting stage 518 based on a 
20 spatial shifting matrix, S, thereby generating shifted upsampled image 520. In 
the illustrated embodiment, shifting stage 518 performs a one pixel diagonal 
shift. Images 504 and 520 are convolved with an interpolating filter at 
convolution stages 506 and 522, respectively, thereby generating blocked 
images. In the illustrated embodiment, the interpolating filter at convolution 
25 stages 506 and 522 is a 2x2 filter with filter coefficients of "1", and with the 
center of the convolution being the upper left position in the 2x2 matrix. The 
blocked images generated at convolution stages 506 and 522 are summed by 
summation block 508, and multiplied by a factor of 0.5 at multiplication stage 
510, to generate the 8x8 pixel simulated high resolution image 512. The image 
30 data is multiplied by a factor of 0.5 at multiplication stage 5 10 because, in one 
embodiment, each of the sub-frames 30F and 30G is displayed for only half of 
the time slot per period allotted to a color. In another embodiment, rather than 
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multiplying by a factor of 0.5 at multiplication stage 510, the filter coefficients 
of the interpolating filter at stages 506 and 522 are reduced by a factor of 0.5. 

In one embodiment, as shown in Figure 8 and described above, the low 
resolution sub-frame data is represented by two separate sub-frames 30F and 
5 30G, which are separately upsampled based on a diagonal sampling matrix (i.e., 
separable upsampling). In another embodiment, as described below with 
reference to Figure 9, the low resolution sub-frame data is represented by a 
single sub-frame, which is upsampled based on a non-diagonal sampling matrix 
(i.e., non-separable upsampling). 

10 Figure 9 is a block diagram illustrating a system 600 for generating a 

simulated high resolution image 610 for two-position processing based on non- 
separable upsampling of an 8x4 pixel low resolution sub-frame 30H according to 
one embodiment of the present invention. System 600 includes quincunx 
upsampling stage 602, convolution stage 606, and multiplication stage 608. 

15 Sub-frame 30H is upsampled by quincunx upsampling stage 602 based on a 
quincunx sampling matrix, Q, thereby generating upsampled image 604. The 
dark pixels in upsampled image 604 represent the thirty-two pixels from sub- 
frame 30H, and the light pixels in upsampled image 604 represent zero values. 
Sub-frame 30H includes pixel data for two 4x4 pixel sub-frames for two-position 

20 processing. The dark pixels in the first, third, fifth, and seventh rows of 

upsampled image 604 represent pixels for a first 4x4 pixel sub-frame, and the 
dark pixels in the second, fourth, sixth, and eighth rows of upsampled image 604 
represent pixels for a second 4x4 pixel sub-frame. 

The upsampled image 604 is convolved with an interpolating filter at 

25 convolution stage 606, thereby generating a blocked image. In the illustrated 
embodiment, the interpolating filter is a 2x2 filter with filter coefficients of "1", 
and with the center of the convolution being the upper left position in the 2x2 
matrix. The blocked image generated by convolution stage 606 is multiplied by 
a factor of 0.5 at multiplication stage 608, to generate the 8x8 pixel simulated 

30 high resolution image 610. 

Figure 10 is a block diagram illustrating a system 700 for generating a 
simulated high resolution image 706 for four-position processing based on sub- 
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frame 301 according to one embodiment of the present invention. In the 
embodiment illustrated in Figure 10, sub-frame 301 is an 8x8 array of pixels. 
Sub-frame 301 includes pixel data for four 4x4 pixel sub-frames for four-position 
processing. Pixels A1-A16 represent pixels for a first 4x4 pixel sub-frame, 
5 pixels B1-B16 represent pixels for a second 4x4 pixel sub-frame, pixels CI -CI 6 
represent pixels for a third 4x4 pixel sub-frame, and pixels D1-D16 represent 
pixels for a fourth 4x4 pixel sub-frame. 

The sub-frame 301 is convolved with an interpolating filter at 
convolution stage 702, thereby generating a blocked image. In the illustrated 

10 embodiment, the interpolating filter is a 2x2 filter with filter coefficients of "1", 
and with the center of the convolution being the upper left position in the 2x2 
matrix. The blocked image generated by convolution stage 702 is multiplied by 
a factor of 0.25 at multiplication stage 704, to generate the 8x8 pixel simulated 
high resolution image 706. The image data is multiplied by a factor of 0.25 at 

15 multiplication stage 704 because, in one embodiment, each of the four sub- 
frames represented by sub-frame 301 is displayed for only one fourth of the time 
slot per period allotted to a color. In another embodiment, rather than 
multiplying by a factor of 0.25 at multiplication stage 704, the filter coefficients 
of the interpolating filter are correspondingly reduced. 

20 

V. Generation of Sub-frames Based on Error Minimization 

As described above, systems 400, 500, 600, and 700 generate simulated 
high resolution images 412, 512, 610, and 706, respectively, based on low 
resolution sub-frames. If the sub-frames are optimal, the simulated high 

25 resolution image will be as close as possible to the original high resolution 
image 28. Various error metrics may be used to determine how close a 
simulated high resolution image is to an original high resolution image, 
including mean square error, weighted mean square error, as well as others. 

Figure 1 1 is a block diagram illustrating the comparison of a simulated 

30 high resolution image 412/5 12/610/706 and a desired high resolution image 28 
according to one embodiment of the present invention. A simulated high 
resolution image 412, 512, 610, or 706, is subtracted on a pixel-by-pixel basis 
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15 



20 



from high resolution image 28 at subtraction stage 802. In one embodiment, the 
resulting error image data is filtered by a human visual system (HVS) weighting 
filter (W) 804. In one form of the invention, HVS weighting filter 804 filters the 
error image data based on characteristics of the human visual system. In one 
embodiment, HVS weighting filter 804 reduces or eliminates low frequency 
errors. The mean squared error of the filtered data is then determined at stage 
806 to provide a measure of how close the simulated high resolution image 412, 
512, 610, or 706 is to the desired high resolution image 28. 

In one embodiment, systems 400, 500, 600, and 700 are represented 
mathematically in an error cost equation that measures the difference between a 
simulated high resolution image 412, 512, 610, or 706, and the original high 
resolution image 28. Optimal sub-frames are identified by solving the error cost 
equation for the sub-frame data that provides the minimum error between the 
simulated high resolution image and the desired high resolution image. In one 
embodiment, globally optimum solutions are obtained in the spatial domain and 
in the frequency domain, and a locally optimum solution is obtained using an 
adaptive multi-pass algorithm. The spatial domain, frequency domain, and 
adaptive multi-pass algorithms are described in further detail below with 
reference to Figures 12-18. 

VI. Spatial Domain 

A spatial domain solution for generating optimal sub-frames according to 
one embodiment is described in the context of the system 600 shown in Figure 9. 
The system 600 shown in Figure 9 can be represented mathematically in an error 
cost function by the following Equation IX: 

Equation IX 



Z* = argmin J = argmin^J £/ Q (k)/(n-k)-/i(n) 




where: 



1*q = optimal low resolution data for sub-frame 



30H; 
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J = error cost function to be minimized; 

n and k = indices for identifying high resolution 

pixel locations for images 604 and 610; 

lo(k) = image data from upsampled image 604 at 
5 location k; 

f(n - k) = filter coefficient of the interpolating 

filter at a position n - k; and 

h(n) = image data for desired high resolution 

image 28 at location n. 
10 The summation of "lo(k)f[n - k)" in Equation EX represents the 

convolution of the upsampled image 604 and the interpolating filter,/, performed 
at stage 606 in system 600. The filter operation is performed by essentially 
sliding the lower right pixel of the 2x2 interpolating filter over each pixel of the 
upsampled image 604. The four pixels of the upsampled image 604 within the 
15 2x2 interpolating filter window are multiplied by the corresponding filter 
coefficient (i.e., "1" in the illustrated embodiment). The results of the four 
multiplications are summed, and the value for the pixel of the upsampled image 
604 corresponding to the lower right position of the interpolating filter is 
replaced by the sum of the four multiplication results. The high resolution data, 
20 h(n), from the high resolution image 28 is subtracted from the convolution value, 
lQ(^)A n - k), to provide an error value. The summation of the squared error 
over all of the high resolution pixel locations provides a measure of the error to 
be minimized. 

An optimal spatial domain solution can be obtained by taking the 
25 derivative of Equation IX with respect to each of the low resolution pixels, and 
setting it equal to zero as shown in the following Equation X: 
Equation X 

dJ =0, te0 



d/g(t) 

where 



0 = the set of quincunx lattice points. 
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Thus, as can be seen from Equation X, the derivative is taken only at the 
set of quincunx lattice points, which correspond to the dark pixels in upsampled 
image 604 in Figure 9. Inserting the equation for J given in Equation IX into 
Equation X, and taking the derivative as specified in Equation X, results in the 
5 following Equation XI: 



The symbol, Cff, in Equation XI represents the auto-correlation 
coefficients of the interpolating filter,/, as defined by the following Equation 
10 XII: 

Equation XII 



Equation XI 



X/G( k ) C #( t - k ) = S^ n )/( n - t )' teQ 



k n 



C # (n) = £/(n)/(n + k) 



k 



Equation XI can be put into vector form as shown in the following 



Equation XIII: 



15 



Equation XIII 




tG0 



where: 



20 



Cff= matrix of auto-correlation coefficients of the 
interpolating filter,/. 

l* Q = vector representing the unknown image data 
for sub-frame 30H, as well as "don't care" data 
(i.e., the image data corresponding to the light 
pixels in upsampled image 604); 
hf = vector representing a filtered version of the 



25 



simulated high resolution image 610 using the 
interpolating filter,/. 
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Deleting the rows and columns corresponding to "don't care" data (i.e., 
the data that is not in the set of quincunx lattice points, 0), results in the 
following Equation XIV: 

Equation XIV 

5 where: 

l£ = vector representing only the unknown image 
data for sub - frame 30H. 

The above Equation XIV is a sparse non-Toeplitz system representing a 
sparse system of linear equations. Since the matrix of auto-correlation 
coefficients is known, and the vector representing the filtered version of the 
10 simulated high resolution image 610 is known, Equation XIV can be solved to 
determine the optimal image data for sub-frame 30H. In one embodiment, sub- 
frame generation unit 36 is configured to solve Equation XIV to generate sub- 
frames 30. 



15 VII. Frequency Domain 

A frequency domain solution for generating optimal sub-frames 30 
according to one embodiment is described in the context of the system 500 
shown in Figure 8. Before describing the frequency domain solution, a few 
properties of the fast fourier transform (FFT) that are applicable to the frequency 

20 domain solution are described with reference to Figures 12 and 13. 

Figure 12 is a diagram illustrating the effect in the frequency domain of 
the upsampling of a 4x4 pixel sub-frame 30J according to one embodiment of 
the present invention. As shown in Figure 12, sub-frame 30J is upsampled by a 
factor of two by upsampling stage 902 to generate an 8x8 pixel upsampled image 

25 904. The dark pixels in upsampled image 904 represent the sixteen pixels from 
sub-frame 30J, and the light pixels in upsampled image 904 represent zero 
values. Taking the FFT of sub-frame 30J results in image (L) 906. Taking the 
FFT of upsampled image 904 results in image (L v ) 908. Image (L v ) 908 
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includes four 4x4 pixel portions, which are image portion (L;) 910A, image 
portion (L 2 ) 910B, image portion (L 3 ) 910C, and image portion (L 4 ) 910D. As 
shown in Figure 12, image portions 910A-910D are each the same as image 906 
(i.e., Lj = L 2 = L 3 = L 4 = L). 
5 Figure 13 is a diagram illustrating the effect in the frequency domain of 

the shifting of an 8x8 pixel upsampled sub-frame 904 according to one 
embodiment of the present invention. As shown in Figure 13, upsampled sub- 
frame 904 is shifted by shifting stage 1002 to generate shifted image 1004. 
Taking the FFT of upsampled sub-frame 904 results in image (Lu) 1006. Taking 

10 the FFT of shifted image 1004 results in image (L V S) 1008. Image (L V S) 1008 
includes four 4x4 pixel portions, which are image portion (LSj) 101 OA, image 
portion (LS 2 ) 101 0B, image portion (LS 3 ) 10 10C, and image portion (LS 4 ) 
1010D. As shown in Figure 13, image 1008 is the same as image 1006 
multiplied by a complex exponential, W, (i.e., L V S = W* Lu) y where "•" denotes 

15 pointwise multiplication. The values for the complex exponential, W 9 are given 
by the following Equation XV: 
Equation XV 

where: 

k\ = row coordinate in the FFT domain; 
20 li2 = column coordinate in the FFT domain; 

M = number of columns in the image; and 
AT = number of rows in the image. 
The system 500 shown in Figure 8 can be represented mathematically in 
an error cost function by the following Equation XVI: 
25 Equation XVI 

{V A , V B ) = ar gmi n J = argmi n £ [f\ (l a + W f L J- H| f [f,^ + W ( L B )-H ( ] 

(L A X B ) (L A X B ) i 

where: 

(L* A , L*b) = vectors representing the optimal FFT's of 
sub-frames 30F and 30G, respectively, shown in Figure 8; 
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/ = error cost function to be minimized; 
i = index identifying FFT blocks that are averaged (e.g., 
for image 908 in Figure 12, four blocks are averaged, 
with z=l corresponding to block 91 OA, i=2 corresponding 
5 to block 91 OB, i=3 corresponding to block 9 10C, and i=4 

corresponding to block 910D); 

F = matrix representing the FFT of the interpolating filter, 

L A = vector representing the FFT of sub-frame 30F shown 
10 in Figure 8; 

L B = vector representing the FFT of sub-frame 30G 
shown in Figure 8; 

W = matrix representing the FFT of the complex 
coefficient given by Equation XV; 
15 H = vector representing the FFT of the desired high 

resolution image 28. 
The superscript "H" in Equation XVI represents the Hermitian (i.e., X H is 
the Hermitian of X). The "hat" over the letters in Equation XVI indicates that 
those letters represent a diagonal matrix, as defined in the following Equation 
20 XVII: 

Equation XVII 



X = diag(X) = 



r X x 0 0 0 ^1 

0 X 2 0 0 

0 0 X 3 0 

0 0 0 xj 



Taking the derivative of Equation XVI with respect to the complex 
conjugate of L A and setting it equal to zero results in the following Equation 
25 XVffl: 
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Equation XVIII 



Taking the derivative of Equation XVI with respect to the complex 
conjugate of L B and setting it equal to zero results in the following Equation 
5 XEX: 

Equation XIX 

= X^&l* + E*A L * - Z wfc = o 



The horizontal bar over the letters in Equations XVIII and XIX indicates 
that those letters represent a complex conjugate (i.e., A represents the complex 
10 conjugate of A). 

Solving Equations XVIII and XIX for L A and L B results in the following 
Equations XX and XXI 

Equation XX 



20 



L s =^A" 1 b)" 1 (d-A- , c) 



15 Equation XXI 

L A =A (C-BLj 

Equations XX and XXI may be implemented in the frequency domain 
using pseudo-inverse filtering. In one embodiment, sub-frame generation unit 
36 is configured to generate sub-frames 30 based on Equations XX and XXL 



Vin. Adaptive Multi-Pass 

An adaptive multi-pass algorithm for generating sub-frames 30 according 
to one embodiment uses past errors to update estimates for the sub-frame data, 
and provides fast convergence and low memory requirements. The adaptive 
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multi-pass solution according to one embodiment is described in the context of 
the system 600 shown in Figure 9. The system 600 shown in Figure 9 can be 
represented mathematically in an error cost function by the following Equation 
XXII: 

5 Equation XXII 

7 w (n) = |^")(n)| 2 = ^/W(k)/(n - k) - h(n) j 
where: 

n = index identifying the current iteration; 
f n) (n) = error cost function at iteration n; 
e M (n) = square root of the error cost function, 
10 f n) (n); 

n and k = indices for identifying high resolution 
pixel locations in images 604 and 610; 
lQ (n) (k) - image data from upsampled image 604 at 
location k; 

15 f(n - k) = filter coefficient of the interpolating 

filter at a position n - k; and 
h(n) = image data for desired high resolution 
image 28 at location n. 
As can be seen from Equation XXII, rather than minimizing a global 
20 spatial domain error by summing over the entire high resolution image as shown 
in Equation IX above, a local spatial domain error, which is a function of n, is 
being minimized. 

A least mean squares (LMS) algorithm is used in one embodiment to 
determine the update, which is represented in the following Equation XXIII: 
25 Equation XXIII 



Cm-^ct)**^. tee 



where: 

0 = the set of quincunx lattice points (i.e., the dark 
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pixels in upsampled image 604 in Figure 9); and 
a = sharpening factor. 
Taking the derivative of Equation XXII provides the value for the 
derivative in Equation XXIII, which is given in the following Equation XXIV: 



In one embodiment, a block-LMS algorithm using the average gradient 
over a "region of influence" is used to perform the update, as represented by the 
following Equation.XXV: 



where: 

Q = region of influence 
Figure 14 is a diagram illustrating regions of influence (Q) 1 106 and 
1 108 for pixels in an upsampled image 1 100 according to one embodiment of the 

15 present invention. Pixel 1 102 of image 1 100 corresponds to a pixel for a first 
sub-frame, and pixel 1 104 of image 1 100 corresponds to a pixel for a second 
sub-frame. Region 1 106, which includes a 2x2 array of pixels with pixel 1 102 in 
the upper left corner of the 2x2 array, is the region of influence for pixel 1 102. 
Similarly, region 1 108, which includes a 2x2 array of pixels with pixel 1 104 in 

20 the upper left corner of the 2x2 array, is the region of influence for pixel 1 104. 

Figure 15 is a diagram illustrating the generation of an initial simulated 
high resolution image 1208 based on an adaptive multi-pass algorithm according 
to one embodiment of the present invention. An initial set of low resolution sub- 
frames 30K-1 and 30L-1 are generated based on an original high resolution 

25 image 28. In the illustrated embodiment, the initial set of sub-frames 30K-1 and 
30L-1 are generated using an embodiment of the nearest neighbor algorithm 
described above with reference to Figure 5. The sub-frames 30K-1 and 30L-1 



Equation XXIV 




Equation XXV 



dJ {n) (n) 
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are upsampled to generate upsampled image 1202, The upsampled image 1202 
is convolved with an interpolating filter 1204, thereby generating a blocked 
image, which is then multiplied by a factor of 0.5 to generate simulated high 
resolution image 1208. In the illustrated embodiment, the interpolating filter 
5 1204 is a 2x2 filter with filter coefficients of "1", and with the center of the 
convolution being the upper left position in the 2x2 matrix. The lower right 
pixel 1206 of the interpolating filter 1204 is positioned over each pixel in image 
1202 to determine the blocked value for that pixel position. As shown in Figure 
15, the lower right pixel 1206 of the interpolating filter 1204 is positioned over 

10 the pixel in the third row and fourth column of image 1202, which has a value of 
"0". The blocked value for that pixel position is determined by multiplying the 
filter coefficients by the pixel values within the window of the filter 1204, and 
adding the results. Out-of- frame values are considered to be "0". For the 
illustrated embodiment, the blocked value for the pixel in the third row and 

15 fourth column of image 1202 is given by the following Equation XXVI 
Equation XXVI 

(1 x 0) + (1 x 5) + (1 x 5) + (1 x 0) = 10 

The value in Equation XXVI is then multiplied by the factor 0.5, and the 
result (i.e., 5) is the pixel value for the pixel 1210 in the third row and the fourth 

20 column of the initial simulated high resolution image 1208. 

After the initial simulated high resolution image 1208 is generated, 
correction data is generated. Figure 16 is a diagram illustrating the generation of 
correction data based on the adaptive multi-pass algorithm according to one 
embodiment of the present invention. As shown in Figure 16, the initial 

25 simulated high resolution image 1208 is subtracted from the original high 

resolution image 28 to generate an error image 1302. Correction sub-frames 
1312 and 1314 are generated by averaging 2x2 blocks of pixels in error image 
1302. For example, the pixel 1308 in the first column and first row of error 
image 1302 has a region of influence 1304. The pixel values within the region 

30 of influence 1304 are averaged to generate a first correction value (i.e., 0.75). 
The first correction value is used for the pixel in the first column and the first 
row of correction sub-frame 1312. Similarly, the pixel 1310 in the second 
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column and second row of error image 1302 has a region of influence 1306. The 
pixel values within the region of influence 1306 are averaged to generate a 
second correction value (i.e., 0.75). The second correction value is used for the 
pixel in the first column and the first row of correction sub-frame 1314. 
5 The correction value in the first row and second column of correction 

sub-frame 1312 (i.e., 1.38) is generated by essentially sliding the illustrated 
region of influence box 1304 two columns to the right and averaging those four 
pixels within the box 1304. The correction value in the second row and first 
column of correction sub-frame 1312 (i.e., 0.50) is generated by essentially 

10 sliding the illustrated region of influence box 1304 two rows down and 

averaging those four pixels within the box 1304. The correction value in the 
second row and second column of correction sub-frame 1312 (i.e., 0.75) is 
generated by essentially sliding the illustrated region of influence box 1304 two 
columns to the right and two rows down and averaging those four pixels within 

15 the box 1304. 

The correction value in the first row and second column of correction 
sub-frame 1314 (i.e., 0.00) is generated by essentially sliding the illustrated 
region of influence box 1306 two columns to the right and averaging those pixels 
within the box 1306. Out-of-frame values are considered to be "0". The 

20 correction value in the second row and first column of correction sub-frame 
1314 (i.e., 0.38) is generated by essentially sliding the illustrated region of 
influence box 1306 two rows down and averaging those pixels within the box 
1306. The correction value in the second row and second column of correction 
sub-frame 1314 (i.e., 0.00) is generated by essentially sliding the illustrated 

25 region of influence box 1306 two columns to the right and two rows down and 
averaging those four pixels within the box 1306. 

The correction sub-frames 1312 and 1314 are used to generate updated 
sub-frames. Figure 17 is a diagram illustrating the generation of updated sub- 
frames 30K-2 and 30L-2 based on the adaptive multi-pass algorithm according 

30 to one embodiment of the present invention. As shown in Figure 17, the updated 
sub-frame 30K-2 is generated by multiplying the correction sub-frame 1312 by 
the sharpening factor, a, and adding the initial sub-frame 30K-1. The updated 
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sub-frame 30L-2 is generated by multiplying the correction sub-frame 1314 by 
the sharpening factor, a, and adding the initial sub-frame 30L-1. In the 
illustrated embodiment, the sharpening factor, a, is equal to 0.8. 

In one embodiment, updated sub-frames 30K-2 and 30L-2 are used in the 
5 next iteration of the adaptive multi-pass algorithm to generate further updated 
sub-frames. Any desired number of iterations may be performed. After a 
number of iterations, the values for the sub-frames generated using the adaptive 
multi-pass algorithm converge to optimal values. In one embodiment, sub-frame 
generation unit 36 is configured to generate sub-frames 30 based on the adaptive 
10 multi-pass algorithm. 

The embodiment of the adaptive multi-pass algorithm described above 
with reference to Figures 15-17 is for two-position processing. For four-position 
processing, Equation XXIV becomes the following Equation XXVII: 

Equation XXVII 



= 2 X^ (rt) (k)/(n-k)-/z(n) 



/(n-t) 



3/ (w) (t) ^ 
15 where: 

f n) = low resolution data for the four sub-frames 
30; 

And Equation XXH[ becomes the following Equation XXVIQ: 
20 Equation XXVIII 

/M(t) = /W(t) + a^^ 

a/ (n) (t) 

For four-position processing, there are four sub-frames, so the amount of 
low resolution data is the same as the amount of high resolution data. Each high 
resolution grid point contributes one error, and there is no need to average 
gradient update as represented in Equation XXV above. Rather, the error at a 
25 given location directly gives the update. 

As described above, in one embodiment, the adaptive multi-pass 
algorithm uses a least mean squares (LMS) technique to generate correction 
data. In another embodiment, the adaptive multi-pass algorithm uses a 
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projection on a convex set (POCS) technique to generate correction data. The 
adaptive multi-pass solution based on the POCS technique according to one 
embodiment is described in the context of the system 600 shown in Figure 9. 
The system 600 shown in Figure 9 can be represented mathematically in an error 
cost function by the following Equation XXIX: 
Equation XXIX 



e(n)\ = 



|J]Z G (k)/(n-k)--/i(n) 
where: 

e(n) = error cost function; 

n and k = indices identifying high resolution pixel 
10 locations; 

lQ(k) = image data from upsampled image 604 at 
location k; 

f(n - k) = filter coefficient of the interpolating 
filter at a position n - k; and 
15 h(n) = image data for desired high resolution 

image 28 at location n. 
A constrained set for the POCS technique is defined by the following 
Equation XXX: 

Equation XXX 

C(n) = |/ e (n): ^l Q (k)f(n-k)-h(n)j <//j 
20 where: 

C(n) = constrained set that includes all sub-frame 
data from upsampled image 604 that is bounded by 
parameter, r\; and 

r| = error magnitude bound constraint. 
25 The sub-frame pixel values for the current iteration are determined based 

on the following Equation XXXI: 
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Equation XXXI 

(l-^)/y(t) + A cv jjy~ f# e(n*)>^ (te0) 

(l-A)/W(t) + ^^f^7 
l/l 

Z< n) (t) e(n) = /7 



/r ) (t) = <! 



where: 

n = index identifying the current iteration; 
5 X = relaxation parameter; and 

||/|| = norm of the coefficients of the interpolating 
filter. 

The symbol, n*, in Equation XXXI represents the location in the region 
of influence, Q, where the error is a maximum, and is defined by the following 
10 Equation XXXII: 

Equation XXXII 

n* = argmaxjn e Q : |e(n)|} 

Figure 18 is a diagram illustrating the generation of correction data based 
on the adaptive multi-pass algorithm using a POCS technique according to one 

15 embodiment of the present invention. In one embodiment, an initial simulated 
high resolution image 1208 is generated in the same manner as described above 
with reference to Figure 15, and the initial simulated high resolution image 1208 
is subtracted from the original high resolution image 28 to generate an error 
image 1302. The Equation XXXI above is then used to generate updated sub- 

20 frames 30K-3 and 30L-3 from the data in error image 1302. For the illustrated 
embodiment, it is assumed that relaxation parameter, A, in Equation XXXI is 
equal to 0.5, and the error magnitude bound constraint, r\ 9 is equal to 1. 

With the POCS technique, rather than averaging the pixel values within 
the region of influence to determine a correction value as described above with 

25 reference to Figure 16, the maximum error, e(n*), within the region of influence 
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is identified. An updated pixel value is then generated using the appropriate 
formula from Equation XXXI, which will depend on whether the maximum 
error, e(n*), within the region of influence is greater than 1, less than 1, or equal 
to 1 (since r| = 1 for this example). 
5 For example, the pixel in the first column and first row of error image 

1302 has a region of influence 1304. The maximum error within this region of 
influence 1304 is 1 (i.e., e(n*) = 1). Referring to Equation XXXI, for the case 
where e(n*) = 1 , the updated pixel value is equal to the previous value for this 
pixel. Referring to Figure 15, the previous value for the pixel in the first column 

10 and the first row of sub-frame 30K-1 was 2, so this pixel remains with a value of 
2 in updated sub-frame 30K-3. The pixel in the second column and second row 
of error image 1302 has a region of influence 1306. The maximum error within 
this region of influence 1306 is 1.5 (i.e., e(n*) = 1.5). Referring to Equation 
XXXI, for the case where e(n*) > 1, the updated pixel value is equal to half the 

15 previous value for this pixel, plus half of the quantity (e(n*) - 1), which is equal 
to 1.25. Referring to Figure 15, the previous value for the pixel in the first 
column and the first row of sub-frame 30L-1 was 2, so the updated value for this 
pixel is 1.25 in updated sub-frame 30L-3. 

The region of influence boxes 1302 and 1304 are essentially moved 

20 around the error image 1302 in the same manner as described above with 

reference to Figure 1 6 to generate the remaining updated values in updated sub- 
frames 30K-3 and 30L-3 based on Equation XXXI. 

IX. Format Conversion 

25 In one embodiment, analog image data 162 (Figure 1) is sampled by A/D 

converter 32 (Figure 1) on a rectangular grid. In the embodiments described 
above, the desired high resolution image 28 and the generated sub-frames 30 are 
made up of rectangular-shaped pixels arranged on rectangular grids. Figure 19A 
is a diagram illustrating rectangular-shaped (e.g., square) pixels 1360 on a 

30 rectangular grid 1362 according to one embodiment of the present invention. 

Four neighboring rectangular-shaped pixels 1360 are shown in Figure 19A. The 
centers 1364 of the pixels 1360 define a rectangular grid 1362. It will be 
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understood by persons of ordinary skill in the art that rectangular grid 1362, 
which is shown with four grid points 1364 and four pixels 1360, may include 
any desired number of grid points 1364 and pixels 1360. 

Also shown in Figure 19A is a pair of orthogonal axes 1366 and 1368. 
5 X-axis 1368 represents a horizontal dimension, and Y-axis 1366 represents a 
vertical dimension. A row of pixels 1360 on rectangular grid 1362 is defined by 
drawing a line through the centers of pixels 1360 parallel to the horizontal 
dimension represented by X-axis 1368. A column of pixels 1360 on rectangular 
grid 1362 is defined by drawing a line through the centers of pixels 1360 parallel 

10 to the vertical dimension represented by Y-axis 1366. The rectangular grid 1362 
shown in Figure 19A includes two rows and two columns of pixels 1360. 

When displayed, each row of pixels 1360 on rectangular grid 1362 is 
parallel to the horizontal dimension represented by X-axis 1368, and each 
column of pixels 1360 on rectangular grid 1362 is parallel to the vertical 

15 dimension represented by Y-axis 1366. Also, any two adjacent pixels 1360 that 
are in the same row will be in adjacent columns, and any two adjacent pixels 
1360 that are in the same column will be in adjacent rows. 

In another embodiment of the present invention, analog image data 162 is 
sampled by A/D converter 32 on a diamond grid. In one embodiment, the 

20 desired high resolution image 28 and the generated sub-frames 30 are made up 
of diamond-shaped pixels arranged on diamond grids. Figure 19B is a diagram 
illustrating diamond-shaped pixels 1370 on a diamond grid 1372 according to 
one embodiment of the present invention. Four neighboring diamond-shaped 
pixels 1370 are shown in Figure 19B. The centers 1374 of the pixels 1370 

25 define a diamond grid 1372. Diamond grid 1372 is also referred to as a 

quincunx grid. It will be understood by persons of ordinary skill in the art that 
diamond grid 1372, which is shown with four grid points 1374 and four pixels 
1370, may include any desired number of grid points 1374 and pixels 1370. 

A row of pixels 1370 on diamond grid 1372 is defined by drawing a line 

30 through the centers of pixels 1370 parallel to the horizontal dimension 

represented by X-axis 1368. A column of pixels 1370 on diamond grid 1372 is 
defined by drawing a line through the centers of pixels 1370 parallel to the 
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vertical dimension represented by Y-axis 1366. The diamond grid 1372 shown 
in Figure 19B includes three rows and three columns of pixels 1370. 

When displayed, each row of pixels 1370 on diamond grid 1372 is 
parallel to the horizontal dimension represented by X-axis 1 368, and each 
5 column of pixels 1370 on diamond grid 1372 is parallel to the vertical dimension 
represented by Y-axis 1366. However, unlike rectangular grid 1362, when 
pixels 1370 on diamond grid 1372 are displayed, any two adjacent pixels 1370 
that are in the same row will not be in adjacent columns, and any two adjacent 
pixels 1370 that are in the same column will not be in adjacent rows. The pixels 

10 1370 in adjacent rows are offset from one another in the horizontal dimension 
represented by X-axis 1368, rather than being aligned with one another like 
pixels 1360 on rectangular grid 1362. 

In one form of the invention, sub-frame generation unit 36 (Figure 1) is 
configured to generate sub-frames 30 having rectangular grids of rectangular- 

15 shaped pixels, and display device 26 (Figure 1) is a low resolution rectangular 
display that is configured to display these generated rectangular sub-frames 30. 
In another form of the invention, sub-frame generation unit 36 is configured to 
generate sub-frames 30 having diamond grids of diamond-shaped pixels, and 
display device 26 is a low resolution quincunx display that is configured to 

20 display these generated diamond sub-frames 30. Diamond grids have some 

advantages over rectangular grids, including: (1) The sub-frame pixels line up in 
rows and columns so it is possible to compute a resolution specification; and (2) 
signals sampled on diamond grids alias first along the diagonal frequencies 
where humans are less visually sensitive, rather than at the vertical and 

25 horizontal frequencies. 

In one embodiment, sub-frame generation unit 36 is configured to 
perform a sampling format to display format conversion function (e.g., 
converting from rectangular grids to diamond grids, or from diamond grids to 
rectangular grids). In one form of the invention, sub-frame generation unit 36 is 

30 configured to generate sub-frames 30 having diamond-shaped pixels on a 

diamond grid based on a high-resolution image 28 sampled on a rectangular grid 
and having rectangular-shaped pixels, wherein the generated sub-frames are 
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displayed using diamond two-position processing. One embodiment of a 
method performed by sub-frame generation unit 36 for converting from a 
rectangular sampling format to a diamond display format is described below 
with reference to Figures 20-22. 
5 In another form of the invention, sub-frame generation unit 36 is 

configured to generate an image having diamond-shaped pixels on a diamond 
grid based on a high-resolution image 28 sampled on a rectangular grid and 
having rectangular-shaped pixels, wherein the generated image is displayed 
without using two-position or four-position processing. One embodiment of a 
10 method performed by sub-frame generation unit 36 for converting from a 

rectangular sampling format to a diamond display format without using two- 
position or four-position processing is described below with reference to Figure 
23. 

In yet another form of the invention, sub-frame generation unit 36 is 

15 configured to generate sub-frames 30 having rectangular- shaped pixels on a 

rectangular grid based on a high-resolution image 28 sampled on a diamond grid 
and having diamond-shaped pixels. One embodiment of a method performed by 
sub-frame generation unit 36 for converting from a diamond sampling format to 
a rectangular display format is described below with reference to Figures 24-25. 

20 Figure 20 is a diagram illustrating the display of two sub-frames 30M 

and 30N with rectangular-shaped pixels on rectangular grids and a diagonal 
offset between sub-frames according to one embodiment of the present 
invention. Sub-frame 30M includes four low resolution rectangular-shaped 
pixels 1402A-1402D (shown in Figure 20 with fine stipple shading) organized in 

25 two rows and two columns, and sub-frame 30N includes four low resolution 

rectangular-shaped pixels 1404A-1404D (shown in Figure 20 with coarse stipple 
shading) organized in two rows and two columns. The pixels for sub-frames 
30M and 30N are each arranged on a low resolution rectangular grid. 

Sub-frame 30N is shifted in a diagonal direction (e.g., down and to the 

30 right) with respect to sub-frame 30M by a sub-pixel amount (e.g., one half 
pixel). When sub-frames 30M and 30N are displayed in relatively quick 
succession using two-position processing, the displayed image appears to the 
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human visual system to be essentially a superposition of the two sub-frames 
30M and 30N having a higher resolution than either of the individual sub-frames 
30M and 30N. The displayed image appears to have high resolution pixels 
1406. The high resolution pixels 1406 are also rectangular-shaped, and are 
5 positioned on a high resolution rectangular grid. Each of the high resolution 

pixels 1406 has an area, A. Each of the low-resolution sub-frame pixels 1402A- 
1402D and 1404A-1404D has an area that is four times the area of the high 
resolution pixels 1406 (i.e., 4 x A). 

Optimal pixel values for sub-frames 30M and 30N can be generated 

10 using the adaptive multi-pass algorithm, which is described above with reference 
to the system 600 shown in Figure 9. For two-position processing, in the case of 
a desired high resolution image 28 on a rectangular grid and sub-frames 30 on a 
rectangular grid, the adaptive multi-pass algorithm according to one embodiment 
uses an averaging over a region of influence as described above, since there are 

15 more pixels in the error image (e.g., error image 1302 shown in Figure 16) than 
in the sub-frames (e.g., sub-frames 30K-1 and 30L-1 shown in Figure 15). The 
error is minimized on every high resolution grid point, so the error image has as 
many pixels as the desired high resolution image 28. Each of the two sub-frames 
30 for two-position processing has one-fourth as many pixels as the high 

20 resolution image 28. Thus, the two sub-frames 30 combined will have half as 
many pixels as the high resolution image 28. Hence, each sub-frame pixel is 
updated by averaging some error values, which can limit image quality. 

In one embodiment, the averaging over a region of influence step in the 
adaptive multi-pass algorithm for two-position processing is eliminated by 

25 generating sub-frames 30 on a diamond grid based on a desired high resolution 
image 28 on a rectangular grid, as described in further detail below with 
reference to Figures 21 and 22. 

Figure 21 is a diagram illustrating the display of two sub-frames 30P and 
30Q with diamond-shaped pixels on diamond grids and a horizontal offset 

30 between sub-frames according to one embodiment of the present invention. 
Sub-frame 30P includes four low resolution diamond-shaped pixels 1452A- 
1452D (shown in Figure 21 with fine stipple shading), and sub-frame 30Q 
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includes four low resolution diamond-shaped pixels 1454A-1454D (shown in 
Figure 21 with coarse stipple shading). The pixels for sub-frames 30P and 30Q 
are each arranged on a low resolution diamond grid. 

Sub-frame 30Q is shifted horizontally to the right with respect to sub- 
5 frame 30P by a sub-pixel amount (e.g., one half pixel). When sub-frames 30P 
and 30Q are displayed in relatively quick succession using two-position 
processing, the displayed image appears to the human visual system to be 
essentially a superposition of the two sub-frames 30P and 30Q having higher 
resolution than either of the individual sub-frames 30P and 30Q. The displayed 
10 image appears to have high resolution pixels 1456. The high resolution pixels 
1456 are also diamond-shaped, and are positioned on a high resolution diamond 
grid. The display of two temporally and spatially shifted sub-frames 30 with 
diamond-shaped pixels on diamond grids is referred to herein as diamond two- 
position processing. 

15 In one embodiment, optimal pixel values for sub-frames 30P and 30Q are 

generated based on the adaptive multi-pass algorithm described above, but 
without averaging over a region of influence. In one form of the invention, 
optimal pixel values are generated by assuming that the high resolution image 
formed by the superposition of two sub-frames with diamond-shaped pixels on 

20 diamond grids (e.g., sub-frames 30P and 30Q) has pixels that are larger than 
those assumed for a high resolution image formed by the superposition of the 
two sub-frames 30M and 30N shown in Figure 20. This assumption results in 
there being as many error image pixels as there are sub-frame pixels. As a 
consequence, the sub-frame pixels can be updated without averaging, resulting 

25 in better image quality. 

As shown in Figure 20, the high resolution pixels 1406 produced by two- 
position processing of sub-frames 30M and 30N were assumed to each have an 
area, A. In contrast, for the sub-frames 30P and 30Q shown in Figure 21, it is 
assumed that these sub-frames produce high-resolution pixels 1458 A and 1458B, 

30 which are on a "hypothetical" rectangular grid formed from the superposition of 
two diamond grids. The rectangular grid is referred to as "hypothetical" because 
a rectangular grid is not actually used for the display, but rather it is assumed to 
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be present for the purpose of calculating optimal sub-frame pixel values. 
Assuming that each of the pixels 1456 has an area, A, then each of the high 
resolution pixels 1458A and 1458B will have an area that is twice as large (i.e., 2 
x A) as the pixels 1456. And each of the high resolution pixels 1458A and 
5 1458B will have an area that is one-half as large as the area of the sub-frame 

pixels 1452A-1452D and 1454A-1454D. In contrast, each of the high resolution 
pixels 1406 (Figure 20) has an area that is one-fourth as large as the area of the 
sub-frame pixels 1402A-1402D and 1404A-1404D. 

As shown in Figure 21, several low-resolution pixels in sub-frames 30P 

10 and 30Q contribute to forming each high resolution pixel 1458A and 1458B. 
For example, pixel 1458A is formed from one-eighth of each of the four pixels 
1452A-1452D in sub-frame 30P, and one-half of the pixel 1454B in sub-frame 
30Q. Similarly, pixel 1458B is formed from one-eighth of each of the four 
pixels 1454A-1454D in sub-frame 30Q, and one-half of the pixel 1452C in sub- 

15 frame 30P. Knowing these contributions is useful in formulating a system for 
generating a simulated high resolution image based on diamond two-position 
processing. In one form of the invention, optimal pixel values for sub-frames 
30P and 30Q are generated based on minimization of an error metric between a 
desired high resolution image 28 (Figure 1 1) on a rectangular grid and a 

20 simulated high resolution image. 

Figure 22 is a block diagram illustrating a system 1500 for generating a 
simulated high resolution image 1508 on a rectangular grid based on sub-frames 
on diamond grids according to one embodiment of the present invention. In the 
embodiment illustrated in Figure 22, sub-frame 30R is an 8x8 array of pixels. 

25 Sub-frame 30R includes pixel data for two 32-pixel sub-frames with diamond 
shaped-pixels on diamond grids. Pixels 1502A represent pixels for a first 32- 
pixel sub-frame on a diamond grid, and pixels 1502B represent pixels for a 
second 32-pixel sub-frame on a diamond grid. Even though the pixels 1502A 
and 1502B represent diamond-shaped pixels for two sub-frames on diamond 

30 grids, such as shown in Figure 21, the pixels 1502A and 1502B are shown in 
sub-frame 30R on a rectangular grid, which may be considered a hypothetical 
rectangular grid as it is used for purposes of calculating the simulated high 
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resolution image 1508, and is not used for display. Similarly, simulated high 
resolution image 1508 includes rectangular-shaped pixels 1510 that are on a 
rectangular grid, but this grid may be considered a hypothetical rectangular grid, 
as the grid is assumed to be formed from the superposition of two diamond grids 
5 as described above with reference to Figure 21. 

The sub-frame 30R is convolved with an interpolating filter at 
convolution stage 1506, thereby generating an 8x8 pixel simulated high 
resolution image 1508. In the illustrated embodiment, the interpolating filter has 
four filter coefficients of "1/8", and one filter coefficient of "1/2". The center of 

10 convolution is the middle position in the filter (i.e., the position with filter 

coefficient of "1/2"). As described above with reference to Figure 21, several 
low-resolution pixels in sub-frames 30P and 30Q contribute to forming each 
high resolution pixel 1458A and 1458B. The layout of the interpolating filter 
and the values for the filter coefficients are chosen based on these contributions 

15 shown in Figure 21 and described above. 

In one embodiment, system 1500 is represented mathematically in an 
error cost equation that measures the difference between the simulated high 
resolution image 1508, and an original high resolution image 28. In one 
embodiment, the simulated high resolution image 1508 and the original high 

20 resolution image 28 are both on rectangular grids. Optimal sub-frames are 
identified by solving the error cost equation for the sub-frame data 30R that 
provides the minimum error between the simulated high resolution image 1508 
and the desired high resolution image 28. In one embodiment, a globally 
optimum solution is obtained in the spatial domain based on the spatial domain 

25 algorithm described above. In another embodiment, a globally optimum solution 
is obtained in the frequency domain based on the frequency domain algorithm 
described above. In yet another embodiment, locally optimum solutions are 
obtained based on the adaptive multi-pass algorithm described above. 

The optimal sub-frame data obtained by solving the error cost equation 

30 will be on a rectangular grid, as indicated by sub-frame 30R in Figure 22. 

However, the pixel data in sub-frame 30R is actually used to generate two sub- 
frames on diamond grids with diamond-shaped pixels, which are displayed by 
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display device 26 using diamond two-position processing. Each rectangular- 
shaped pixel 1502A corresponds to a diamond-shaped pixel in a first sub-frame, 
and each rectangular-shaped pixel 1502B corresponds to a diamond-shaped pixel 
in a second sub-frame. Each row of rectangular-shaped pixels 1502A 
5 corresponds to a row of diamond-shaped pixels in the first sub-frame, and each 
row of rectangular-shaped pixels 1502B corresponds to a row of diamond- 
shaped pixels in the second sub-frame. The appearance of the actual display of 
either of the two sub-frames can be shown by drawing diamonds around the 
rectangular-shaped pixels 1502A or 1502B for the sub-frame, as indicated by the 

10 diamonds 1504 drawn with hidden lines around six of the rectangular- shaped 

pixels 1502B. The rectangular-shaped pixels 1502A and 1502B on a rectangular 
grid are used for computing the optimum pixel values. These pixel values are 
then assigned to the diamond-shaped pixels on diamond grids that are actually 
displayed by display device 26. 

15 As shown in Figure 22, there are as many sub-frame pixels 1502A and 

1502B as there are pixels 1510 in simulated high resolution image 1508 (and 
correspondingly in the desired high resolution image 28). Thus, by performing 
the rectangular sampling format to diamond display format conversion, and 
generating optimal sub-frame pixel values based on the adaptive multi-pass 

20 algorithm, the errors do not have to be averaged, which provides a large image 
quality improvement over the display of a rectangular sampled signal on a 
rectangular display using two-position processing of sub-frames on rectangular 
grids. 

As described above with reference to Figures 21 and 22, in one 
25 embodiment, sub-frame generation unit 36 is configured to generate optimal sub- 
frames on a diamond grid for diamond two-position processing based on a high 
resolution image on a rectangular grid. In another embodiment of the present 
invention, sub-frame generation unit 36 is configured to generate an optimal 
image on a diamond grid from a high resolution image on a rectangular grid. In 
30 one embodiment, the generated optimal image on a diamond grid is displayed by 
display device 26 in a conventional manner without using two-position 
processing or four-position processing. In one form of the invention, the pixel 
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values for the optimal image are generated based on minimization of an error 
metric between a desired high resolution image 28 (Figure 1 1) on a rectangular 
grid and a simulated high resolution image. 

Figure 23 is a block diagram illustrating a system 1518 for generating a 
5 simulated high resolution image 1528 on a rectangular grid based on an image 
on a diamond grid according to one embodiment of the present invention. In the 
embodiment illustrated in Figure 23, image 1520 is an 8x8 array of pixels. The 
thirty-two pixels 1522 (shown in Figure 23 with stipple shading) represent pixels 
for an image on a diamond grid, and the remaining thirty-two pixels (shown with 

10 no shading) have a value of zero. Even though the pixels 1522 represent 
diamond-shaped pixels for an image on a diamond grid, the pixels 1522 are 
shown in image 1520 on a rectangular grid, which may be considered a 
hypothetical rectangular grid as it is used for purposes of calculating the 
simulated high resolution image 1528, and is not used for display. 

15 The image 1520 is convolved with an interpolating filter at convolution 

stage 1526, thereby generating an 8x8 pixel simulated high resolution image 
1528. In the illustrated embodiment, the interpolating filter has four filter 
coefficients of "1/2", and one filter coefficient of "1". The center of convolution 
is the middle position in the filter (i.e., the position with filter coefficient of "1"). 

20 In one embodiment, system 1518 is represented mathematically in an 

error cost equation that measures the difference between the simulated high 
resolution image 1528, and an original high resolution image 28. In one 
embodiment, the simulated high resolution image 1528 and the original high 
resolution image 28 are both on rectangular grids. Optimal pixel values for 

25 image 1 520 are identified by solving the error cost equation for the pixel values 
in image 1520 that provide the minimum error between the simulated high 
resolution image 1528 and the desired high resolution image 28. In one 
embodiment, a globally optimum solution is obtained in the spatial domain 
based on the spatial domain algorithm described above. In another embodiment, 

30 a globally optimum solution is obtained in the frequency domain based on the 
frequency domain algorithm described above. In yet another embodiment, 
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locally optimum solutions are obtained based on the adaptive multi-pass 
algorithm described above. 

The optimal pixel data obtained by solving the error cost equation will be 
on a rectangular grid, as indicated by image 1520 in Figure 23. However, the 
5 pixel data in image 1520 is actually used to generate an image on a diamond grid 
with diamond-shaped pixels, which is displayed by display device 26. Each 
rectangular-shaped pixel 1522 corresponds to a diamond-shaped pixel in the 
image to be displayed. Each row of pixels 1522 corresponds to a row of 
diamond-shaped pixels in the image to be displayed. The appearance of the 

10 actual display of the image formed from pixels 1522 can be shown by drawing 
diamonds around each of the rectangular-shaped pixels 1522, as indicated by the 
diamonds 1524 drawn with hidden lines around six of the rectangular-shaped 
pixels 1522. The rectangular- shaped pixels 1522 on a rectangular grid are used 
for computing the optimum pixel values. These pixel values are then assigned to 

15 the diamond-shaped pixels on a diamond grid that are actually displayed by 
display device 26. 

In one embodiment, the averaging over a region of influence step in the 
adaptive multi-pass algorithm for two-position processing is eliminated by 
generating sub-frames 30 on a rectangular grid based on a desired high 

20 resolution image 28 on a diamond grid, as described in further detail below with 
reference to Figures 24 and 25. 

Figure 24 is a diagram illustrating the display of two sub-frames 30S and 
30T with rectangular-shaped pixels on rectangular grids and a diagonal offset 
between sub-frames according to one embodiment of the present invention. 

25 Sub-frame 30S includes four low resolution rectangular-shaped pixels 1552A- 
1552D (shown in Figure 24 with fine stipple shading), and sub-frame 30T 
includes four low resolution rectangular-shaped pixels 1554A-1554D (shown in 
Figure 24 with coarse stipple shading). The pixels for sub-frames 30S and 30T 
are each arranged on a low resolution rectangular grid. 

30 Sub-frame 30T is shifted diagonally (e.g., down and to the right) with 

respect to sub-frame 30S by a sub-pixel amount (e.g., one half pixel). When 
sub-frames 30S and 30T are displayed in relatively quick succession using two- 
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position processing, the displayed image appears to the human visual system to 
be essentially a superposition of the two sub-frames 30S and 30T having higher 
resolution than either of the individual sub-frames 30S and 30T. The displayed 
image appears to have high resolution pixels 1556. The high resolution pixels 
5 1556 are also rectangular-shaped, and are positioned on a high resolution 
rectangular grid. 

In one embodiment, optimal pixel values for sub-frames 30S and 30T are 
generated based on the adaptive multi-pass algorithm described above, but 
without averaging over a region of influence. In one form of the invention, 

10 optimal pixel values are generated by assuming that the high resolution image 
formed by the superposition of two sub-frames with rectangular-shaped pixels 
on rectangular grids (e.g., sub-frames 30S and 30T) has pixels that are larger 
than those assumed for a high resolution image formed by the superposition of 
the two sub-frames 30M and 30N shown in Figure 20. This assumption results 

15 in there being as many error image pixels as there are sub-frame pixels. As a 
consequence, the sub-frame pixels can be updated without averaging, resulting 
in better image quality at a lower computational cost. 

As shown in Figure 20, the high resolution pixels 1406 produced by two- 
position processing of sub-frames 30M and 30N were assumed to each have an 

20 area, A. In contrast, for the sub-frames 30S and 30T shown in Figure 24, it is 

assumed that these sub-frames produce high-resolution pixels 1558A and 1558B, 
which are on a "hypothetical" diamond grid formed from the superposition of 
two rectangular grids. The diamond grid is referred to as "hypothetical" because 
a diamond grid is not actually used for the display, but rather it is assumed to be 

25 present for the purpose of calculating optimal sub-frame pixel values. Assuming 
that each of the pixels 1556 has an area, A, then each of the high resolution 
pixels 1558A and 1558B will have an area that is twice as large (i.e., 2 x A) as 
the pixels 1556. And each of the high resolution pixels 1558A and 1558B will 
have an area that is one-half as large as the area of the sub-frame pixels 1552A- 

30 1552D and 1554A-1554D. In contrast, each of the high resolution pixels 1406 
(Figure 20) has an area that is one-fourth as large as the area of the sub-frame 
pixels 1402A-1402D and 1404A-1404D. 
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As shown in Figure 24, several low-resolution pixels in sub-frames 30S 
and 30T contribute to forming each high resolution pixel 1558A and 1558B. For 
example, pixel 1558A is formed from one-eighth of each of the four pixels 
1552A-1552D in sub-frame 30S, and one-half of the pixel 1554A in sub-frame 
5 30T. Similarly, pixel 1558B is formed from one-eighth of each of the four 

pixels 1554A-1554D in sub-frame 30T, and one-half of the pixel 1552D in sub- 
frame 30S. Knowing these contributions is useful in formulating a system for 
generating a simulated high resolution image based on rectangular two-position 
processing. In one form of the invention, optimal pixel values for sub-frames 

10 30S and 30T are generated based on minimization of an error metric between a 
desired high resolution image 28 (Figure 1 1) on a diamond grid and a simulated 
high resolution image. 

Figure 25 is a block diagram illustrating a system 1600 for generating a 
simulated high resolution image 1608 on a diamond grid based on sub-frames on 

15 rectangular grids according to one embodiment of the present invention. In the 
embodiment illustrated in Figure 25, sub-frame 30U is an 8x8 array of pixels. 
Sub-frame 30U includes pixel data for two 16-pixel sub-frames with rectangular 
shaped-pixels on rectangular grids. Pixels 1602A represent pixels for a first 16- 
pixel sub-frame on a rectangular grid, and pixels 1602B represent pixels for a 

20 second 16-pixel sub-frame on a rectangular grid. The remaining pixels in sub- 
frame 30U, which are indicated by no shading, represent pixels having a value of 
zero. 

The sub-frame 30U is convolved with an interpolating filter at 
convolution stage 1606, thereby generating a simulated high resolution image 

25 1608 with 32 diamond-shaped pixels 1610 on a diamond grid. In the illustrated 
embodiment, the interpolating filter has four filter coefficients of "1/8", and one 
filter coefficient of "1/2". The center of convolution is the middle position in the 
filter (i.e., the position with filter coefficient of "1/2"). As described above with 
reference to Figure 24, several low-resolution pixels in sub-frames 30S and 30T 

30 contribute to forming each high resolution pixel 1558A and 1558B. The layout 
of the interpolating filter and the values for the filter coefficients are chosen 
based on these contributions shown in Figure 24 and described above. 
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In one embodiment, system 1600 is represented mathematically in an 
error cost equation that measures the difference between the simulated high 
resolution image 1608, and an original high resolution image 28. In one 
embodiment, the simulated high resolution image 1608 and the original high 
5 resolution image 28 are both on diamond grids. Optimal sub-frames are 

identified by solving the error cost equation on a diamond grid for the sub-frame 
data 30U that provides the minimum error between the simulated high resolution 
image 1608 and the desired high resolution image 28. In one embodiment, a 
globally optimum solution is obtained in the spatial domain based on the spatial 
10 domain algorithm described above. In another embodiment, a globally optimum 
solution is obtained in the frequency domain based on the frequency domain 
algorithm described above. In yet another embodiment, locally optimum 
solutions are obtained based on the adaptive multi-pass algorithm described 
above. 

15 The optimal sub-frame data obtained by solving the error cost equation is 

used to generate two 16-pixel sub-frames on rectangular grids with rectangular- 
shaped pixels, which are displayed by display device 26 using two-position 
processing. There are as many sub-frame pixels 1602A and 1602B as there are 
pixels 1610 in simulated high resolution image 1508 (and correspondingly in the 

20 desired high resolution image 28). Thus, by performing the diamond sampling 
format to rectangular display format conversion, and generating optimal sub- 
frame pixel values based on the adaptive multi-pass algorithm, the errors do not 
have to be averaged, and the data rate required to produce the same resolution as 
a comparable rectangular sampled image is cut in half, while providing 

25 comparable image quality. Because the data rate is cut in half, less memory is 
needed and a lesser number of computations need to be performed than for a 
comparable rectangular sampled image that is displayed on a rectangular display. 

Although specific embodiments have been illustrated and described 
herein for purposes of description of the preferred embodiment, it will be 

30 appreciated by those of ordinary skill in the art that a wide variety of alternate or 
equivalent implementations may be substituted for the specific embodiments 
shown and described without departing from the scope of the present invention. 
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Those with skill in the mechanical, electro-mechanical, electrical, and computer 
arts will readily appreciate that the present invention may be implemented in a 
very wide variety of embodiments. This application is intended to cover any 
adaptations or variations of the preferred embodiments discussed herein. 
Therefore, it is manifestly intended that this invention be limited only by the 
claims and the equivalents thereof. 
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